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Abstract. We have constructed and modeled a simple and efficient source of slow atoms. From a background 
vapour loaded magneto-optical trap, a thin laser beam extracts a continuous jet of cold rubidium atoms. 
In this setup, the extraction column that is typical to leaking MOT systems is created without any optical 
parts placed inside the vacuum chamber. For detailed analysis, we present a simple 3D numerical simulation 
of the atomic motion in the presence of multiple saturating laser fields combined with an inhomogeneous 
magnetic field. At a pressure of Pjibsr = 1 x 10~*mbar, the moderate laser power of 10 mW per beam 
generates a jet of fiux <P — 1.3 x 10* atoms/s with a mean velocity of 14 m/s and a divergence of < 20 mrad. 



PACS. 32.80.Lg Mechanical effects of light on atoms, molecules, and ions 
atoms; trapping 



32.80.Pj Optical cooling of 



1 Introduction 

Experiments on trapped cold atom clouds require in most 
cases high particle numbers and long trapping lifetimes. In 
order to restrict the lifetime limiting collisions with back- 
ground gas, an ultra-high vacuum (UHV) environment is 
necessary. In turn, at these pressures a purely background 
vapour charged magneto-optical trap (VCMOT) is limited 
to very small atom numbers and long loading times and 
thus needs to be loaded by an additional jet of cold atoms. 

As to the simplest possible cold atom sources, a laser- 
free velocity filter |^ is elegant, but its maximum flux can 
be greatly improved upon by adding a laser cooling stage. 

The Zeeman slower is a widely used technique espe- 
cially for light and thus thermally non capturable fast 
species. For heavier elements, one can accumulate atoms 
into a MOT in a vapour cell, with various strategies for 
subsequent transfer to a recapture MOT in the UHV cell. 
These strategies can be categorized into either a pulsed 
or continuous transfer scheme. The latter category 
involves either a moving molasses or a 'leaking MOT' 
scheme 

This paper presents the construction and numerical 
modeling of a cold atom jet whose flux is continuous, ad- 
justable in a given direction, and velocity tunable. The 
device we present is based on an ordinary VCMOT. It 
captures and cools atoms from the low velocity part of 
the room temperature Maxwell-Boltzmann distribution in 
a high pressure cell of P 10"® mbar. From the center of 
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this source MOT, an additional pushing beam of ~ 1 mm 
spot size extracts a continuous jet that is slow enough 
to be recaptured in a MOT in the UHV region. The jet 
passes through a tube that maintains the pressure differ- 
ential between the two cells, and the atom number transfer 
between the two MOTs is found to be typically 50 % and 
as high as 60 % efficient. 

The Atom Faucet is closely related to the UVIS |^ 
and the 2D+M0T The common concept which relates 
them in the 'leaking MOT' family is the creation of a thin 
extraction column in the center of the MOT where the 
radiation pressure is imbalanced and through which leaks 
a continuous jet of cold atoms. Operation in a continu- 
ous mode maximizes the mean flux up to a value ideally 
equal to the source trap capture rate. Since a leaking trap 
operates at a low trap density, once captured, an atom 
has much higher probability to leave the trap via the jet 
rather than undergoing a collision that expels it. 

The UVIS and 2D+M0T place a mirror inside the vac- 
uum for retroreflection of one of the MOT beams. By 
piercing a hole in this mirror, one creates a hollow retrore- 
flection beam, and the jet exits through the hole. By con- 
trast, the Atom Faucet requires no optical parts inside the 
vacuum system. Here, we superimpose an additional colli- 
mated 'pushing beam' that pierces the extraction column 
through the MOT. 

In these complex magneto-optical arrangements the 
behavior of the system is no longer intuitively obvious. 
On its way into the jet, a thermal atom undergoes sub- 
sequent phases of strong radiation pressure (capture from 
vapour) , overdamped guidance to the magnetic field min- 
imum (MOT molasses) and ID strong radiation pressure 
with transverse 2D molasses cooling (extraction process). 
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Theoretical estimates for near-resonant atom traps con- 
centrate either on the capture Q or on the cooUng 
We develop a simple and heuristic generalization of the 
semiclassical radiation pressure expression for the case of 
multiple saturating laser fields and inhomogeneous mag- 
netic field. The new approach of integrating the atomic 
trajectory through both capture and cooling mechanisms 
(neglecting optical pumping and particle interaction) re- 
produces the parameter dependences of the Atom Faucet. 
The trajectories indicate the physical mechanisms of the 7- 
beam-interplay. However, the simplifications made to the 
Rubidium level scheme lead to an overestimation of the 
absolute value of the radiation pressure force and hence 
an overestimate for the capture velocity of the MOT. 

This paper is organized as follows: In section 2 we 
give details on the experimental realization of the Atom 
Faucet. In section 3 we present the numerical model. Sec- 
tion 4 discusses the parameter dependences of the device 
in the experiment and in the simulations and finally in sec- 
tion 5 we compare this scheme to other vapour cell cold 
atom sources. 



2 Experimental Realisation 

The vacuum system consists of two glass cells separated 
vertically by 67 cm with a MOT aligned at the center of 
each cell. Using an appropriate pumping scheme and a 
differential pumping tube of diameter 5 mm and length 
15 cm the pressure in the lower recapture cell is less than 
10~^^ mbar while in the source cell it is ~ 10~^ mbar. We 
deduce the ®^Rb pressure in the source cell from the reso- 
nant absorption of a multi-passed probe beam. A heated 
reservoir connected to the upper source cell supplies the 
Rubidium vapour. 

A grating stabilized diode laser locked to the |5Si/2, F ~ 
2) |5P3/2,^ — 3) transition injects into three slave 
lasers, two for the source MOT and one for the recapture 
MOT. The Atom Faucet (see fig. 1) is based on a standard 
MOT configuration: two Anti-Hclmholtz-coils maintain a 
magnetic field gradient of 15 G/cm along their axis, which 
is horizontal in this setup. A pair of axial beams with pos- 
itive helicity counterpropagate along the axis of the coils 
and two mutually orthogonal pairs of radial beams with 
negative helicity counterpropagate in the symmetry plane 
of the coils. The radial beams are inclined by 45°. The ra- 
dial trap beams have an 8 mm spot size and the axial beam 
11 mm respectively, all clipped to a diameter of 1 inch by 
our quarterwaveplates. The axial beam carries 20 mW be- 
fore retroreflection, and the radial beams each have 5 mW 
each before retroreflection. The repumping light on the 
|5Si/2,F = 1) ^ |5P3/2,^' = 2) transition from an inde- 
pendent grating stabilized laser is mixed only in the axial 
beam and has a power of 5 mW. 

In addition to these trapping lasers, a permanent push- 
ing beam on the |5Si/2,^' = 2) ^ |5P3/2,F — 3) tran- 
sition with linear polarization and optimal power of 
200 /zW is aligned vertically onto the trap. It is focused to 
a waist of 90 ^m 30 cm before entering the source cell such 
that it diverges to a size of 1.1 mm at the source trap and 




differential pumping tube 




recapture MOT 
in UHV cell 



Fig. 1. The Atom Faucet setup (with the recapture MOT 
below) . A permanent pushing beam with ~ 1 mm spot size 
pierces an extraction column into an ordinary vapour charged 
MOT. The high pressure region is separated from the ultra- 
high- vacuum region by a differential pumping tube. The pres- 
sure in the source cell is monitored by the absorption of an 
additional multi-passed probe beam (not shown). 

3.3 mm at the recapture trap. Its intensity at the center of 
the source MOT and detuning are comparable to those of 
the MOT beams and hence its radiation pressure is also 
comparable with the trapping forces in the MOT. Because 
of the divergence of the pushing beam, the intensity in the 
lower MOT is lower by a factor of 10. It decenters the re- 
capture MOT by ~ 1 mm but does not destabilize it. Note 
that the pushing beam carries no repumping light, so that 
it acts on the atoms only where it intersects the MOT 
beams. 

By studying the loading characteristics of the recap- 
ture MOT, we deduce the main features of the atom jet: 

— When the recapture MOT is empty the initial recap- 
ture rate gives directly the recaptured flux since the 
density dependent intrinsic losses in the MOT are not 
yet important. The absolute number of atoms is deter- 
mined using an absorption imaging technique. 

— The time dependence of the recapture loading rate pro- 
vides a measurement of the longitudinal velocity dis- 
tribution of the jet. More precisely, by suddenly disin- 
jecting the source MOT slave lasers and then recording 
the recapture filling rate via the fluorescence, the char- 
acteristics of the tail of the moving extraction column 



W. Wohlleben, F. Chevy, K. Madison, J. Dalibard,: An Atom Faucet 



3 



> 

O 

(D 
O 

tn 
o 



3.5 



3.0 



2.5 



O 2.0 



CD 



t 1.5 
o 

CD 




20 40 60 80 100 
Time since source deinjection (ms) 



Fig. 2. Development of the fluorescence of the recapture trap 
(circles are photodiode signal) after sudden disinjection of the 
source MOT beams. The pushing beam is not changed in order 
to keep constant its influence on the lower trap fluorescence. 
The fit (solid line) ${v) = x exp ( — (u — w)^/2(5ii^) with a 
Gaussian envelope for the jet velocity distribution yields v — 
14±9m/s. 



are measured. The jet transfer distance D = 67 cm and 
the time delay T of the filling rate response gives the 
mean longitudinal velocity v = D/T in the jet, and 
the time width At of this response gives access to the 
longitudinal velocity dispersion 5v (see Fig. 2). 

For the determination of the transfer efficiency, the 
loading rate of the source MOT is determined by its fluo- 
rescence and compared with the measured recapture flux. 
The fluorescence measurement is done at resonance and 
we assume full saturation of the transition under the in- 
fluence of all six laser beams and thus a photon scattering 
rate of P/2 photons/atom/second. 

We observe a typical transfer efliciency of 50 % (see 
below). Since the radius of the recapture MOT beams is 
r = 5 mm and the transfer distance is Z? = 67 cm, less than 
50% of the atoms are emitted with a divergence larger 
than r/D ^ lOmrad. 



3 Theoretical Description for Numerics 

In order to model both the capture of the atoms from 
the vapor into the source MOT and the subsequent cool- 
ing and pushing processes, we have developed a numeri- 
cal simulation which integrates the equation of motion for 
atoms chosen with random initial positions and velocities. 
We describe the atomic motion using classical dynamics. 
The action of the seven laser beams (6 MOT beams -t- 1 
pushing beam) on an atom located at r with velocity v 
is taken into account through an average radiation force 
F(r, v). We neglect any heating or diffusion caused by 
spontaneous emission. 

The calculation of the semi-classical force acting on an 
atom in this multiple beam configuration is a priori very 



complex. For simplicity, we model the atomic transition as 
a |g, Jg = 0) <-> |e, Je = 1) transition with frequency huiAi 
where \g) and |e) stand for the ground and excited state 
respectively. We denote the lifetime of e. Consider a 
single plane-wave beam with wave vector k, detuning 5 = 
ujl — ujA, intensity J, and polarisation a± along the local 
magnetic field B in r. The radiation pressure force [ p^ 
reads 



2 l-|-s(r,v) 
where the saturation parameter is given by 



(1) 



s(r,v) 



I 



/sat ^2 + 4(5 - k • V ± ^iB/hY ■ 



/i is the magnetic moment associated with level |e) and 
/sat is the saturation intensity for the transition (/^at = 
1.62 mW/cm^ for the D2 resonance line in Rb). Still re- 
stricting our attention to a single traveling wave, we con- 
sider now the case where the light couples \g) to two or 
three Zeeman sublevels |em). The calculation is in this case 
more involved since the solution of the optical Bloch equa- 
tions requires the study of 16 coupled differential equa- 
tions. A simple approximation is obtained in the low sat- 
uration limit [s <^ 1): 



with 



m=-l,0,l 



i(i",v) 



(2) 



/sat r2 + 4(5 - k • V + m/lS/?i)2 

and where Im is the intensity of the laser wave driving 
the \g) jcm) transition. We can sum up the three forces 
associated with the three possible transitions, each calcu- 
lated with the proper detuning taking into account the 
Zeeman effect. 

Still working in the low intensity limit, we can gener- 
alize eq. (^ to the case where N laser beams with wave 
vectors kj and detunings (5j, (j = 1,...,N) are present. 
The force then reads 



.(r,v) 



(3) 



with 



/ 



J, "I 



/sat + 4(<5j - kj • V + m^iB/ny 



Note that in establishing eq. (H) we have taken the 
spatial average of the radiative force over a cell of size 
A = 27r/fc, neglecting thus all interference terms varying 



as i(kj 



r. We therefore neglect any effect of the 



dipole force associated with the light intensity gradients 
on the wavelength scale. This is justified in the case of a 
leaking MOT since the associated dipole potential wells 
are much less deep than the expected residual energy of 
the atoms before extraction. 
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At the center of the capture MOT, we can no longer 
neglect inter-beam saturation effects since the saturation 
parameter for each of the 7 beams is equally ~ 1/7. In 
principle, accounting for this saturation effect requires a 
step- by-step numerical integration of the 16 coupled Bloch 
optical equations (for a \g,Jg = 0) ^ |e, Je = 1) transi- 
tion), as the atom moves in the total electric field result- 
ing from the interference of all the laser beams present in 
the experiment. Such a calculation is unfortunately much 
too computationally intensive to lead to interesting pre- 
dictions for our Atom Faucet in a reasonable time. We 
therefore decided to turn to a heuristic and approximate 
expression for the force, demanding: 

~ In the case of a single traveling wave, a± polarized 
along the magnetic field, we should recover expression 

®- . . . . 

— In the low intensity limit, the force should simplify to 

expression (^. 

— The magnitude of the force should never exceed hk r/2, 
which is the maximal radiation pressure force in a sin- 
gle plane wave. 

There are of course an infinite number of expressions which 
fulfill these three conditions. We have taken the simplest 
one: 

F = ^ ?ik, - (4) 

with partial saturation parameters Sj^m as defined in eq. 
@. This equation is the generalization of the heuristic 
expression used by Phillips and co-workers to account 
for saturation effects in an optical molasses. 

In the simulation, the MOT beams are chosen to have 
Gaussian profiles truncated to the diameter of the quar- 
terwaveplates. Also they are chosen to be equally strong 
with a central intensity of 5 Jgat and to have the proper 
polarizations and directions. The pushing beam's intensity 
is of the same order. We assume that because of optical 
pumping into the lower hyperfine ground state, an atom 
sees no forces when it is out of the repumper light mixed 
in the axial beams. Finally, the magnetic quadrupole field 
is B(x) = b'{-2x,y,z). 

In the simulation the initial position of each atom is 
chosen on one of the cell windows following a uniform spa- 
tial distribution. The initial velocity is given by a Maxwell 
Boltzmann distribution for T = 300 K. The trajectory is 
then integrated using a Runge-Kutta method. From these 
trajectories (see fig. 3), one obtains a probability for an 
atom to be captured and transferred into the jet, as well as 
the jet's characteristics: velocity distribution, divergence, 
and total flux. The absolute flux of the simulated jet is cal- 
ibrated using the real number of atoms emitted per unit 
time and per unit surface of the cell at a pressure P which 
is P/V27rmfcBT |ll|,0. 

The simulation neglects interaction effects like colli- 
sions and multiple light scattering. The validity of the lin- 
ear scaling with pressure is limited to the low pressure 
regime {P < 10~^mbar) where the characteristic extrac- 
tion time of ~ 20 ms is shorter than the collision time, 
which is in turn of the order of the trap lifetime. 




Fig. 3. Some simulated trajectories of atoms in the VCMOT 
+ pushing beam light field that are captured and transfered to 
the jet (distances in mm). 

4 Results 

Inspecting qualitatively the trajectories, we find that an 
atom that enters the beam intersection is first decelerated 
by radiation pressure on a distance much smaller than the 
trapping beam radius. It then slowly moves to the center 
of the trap where it enters the extraction column. The final 
transverse cooling of the jet takes place during extraction, 
so that the divergence of the jet grows if the extraction 
happens too fast. We believe that this is the principal loss 
mechanism of any leaking MOT system, which have in 
common an extraction column and a transverse molasses 
provided by the trapping beams. 



4.1 Total Flux 

For a typical choice of parameters, the simulation finds 
90 % transfer from the source MOT through the differen- 
tial pumping tube to the recapture MOT. The remaining 
10 % of the atoms leave the source at a divergence too large 
to be recaptured and are lost. Experimentally, we have 
achieved a transfer efficiency of at most 60 ± 10%. This 
value is most probably limited by the differential pumping 
tube diameter. 
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Fig. 4. Recaptured flux versus source ceU pressure. The hnear 
fit yields 'P'^^p = 1.3 ± 0.2 x 10** atoms/s x PRb87(10"* mbar). 



then retroreflected. We calculate an optimal detuning of 
sr. The capture rate is divided by more than 2 when the 
detuning is beyond —4.5 r or smaller than —1.5 r. This is 
the typical MOT operation range. The magnetic gradient 
seems to have little influence as long as it is between 8 
and 20G/cm. 

It is particularly helpful to calculate the optimal beam 
waist for a given laser power since in the optical setup this 
parameter is tiresome to change and demands subsequent 
trap realignment. In our case, a 9 mm spot size gives the 
best simulated capture rate, with half maximum values at 
4 mm and 16 mm. For a fixed laser power, having a large 
intersection volume is preferable to increasing the satura- 
tion beyond ~ 4 Igat ■ The experiment uses an 8 mm spot 
size, and the optimum parameters do not change signifi- 
cantly if the retroreflection loss of 20 % is included. Finally, 
the simulation reproduces the smoothly decreasing slope 
of the capture rate versus the MOT beam power of ref M . 



Concerning the total flux, we explored the pressure 
regime of 10~^ < P < 4 x 10~^ mbar and found no devia- 
tion from a linear dependence (see fig. 4) 

<^°^p = 1.3 ± 0.4 X 10* atoms/s x PRb87(10"* mbar). 

The uncertainty primarily comes from the atom number 
determination in the recapture MOT by absorption imag- 
ing. Deviation from linear scaling with pressure is to be 
expected when the collision time with background gas be- 
comes of the order of the typical extraction time from the 
MOT center into the differential pumping tube. This will 
be the case for PnhST > 10^^ mbar. 

In comparison we found that the simulation overesti- 
mates the capture velocity of the MOT, so that we need 
to calibrate its predictions. Therefore we simulate a pure 
MOT without pushing beam and compare the predicted 
capture rate of 

T^^T ^ X IQS atoms/s x PRb87(10"* mbar) 

with the value we measured in the initial regime of linear 
growth of the vapour charged source MOT, 

r^^T ^ 2.5 ± 0.6 X 10* atoms/s x PRb87(10"* mbar) 

We believe that the disagreement between these two 
results corresponds to an overestimation of the source 
MOT capture velocity Vc- Since the number of atoms cap- 
tured in a VCMOT varies as v^, our simple model over- 
estimates Vc by (13/2.5)^/'' ~ 1.5. In the graphs 5,6,7, we 
normalize the absolute value of the flux and concentrate 
on its variation with system parameters. 

Simulated VCMOT Optimisation. Using the simulation 
of a pure MOT without pushing beam, we can readily find 
the parameters which optimise the capture rate from the 
background vapour. The total laser power is taken to be 
20 mW, equally distributed among three beams which are 



4.2 Pushing Beam Parameters 

We now add the thin pushing beam to the MOT light field. 
Doing so does not modify the optimal parameters of the 
capture MOT, neither in experiment nor in the simulation. 
Remember that the volume affected by the thin beam is 
very small compared to the total capture volume of the 
source MOT. We investigate the influence of pushing beam 
power, detuning, and size on the atomic jet emerging from 
the MOT. The following discussion shall directly combine 
experimental findings and the results from the theoretical 
model. 



Power. For very low pushing beam power the trap is de- 
centered but not yet leaking. At Ppush = 80 fiW (cor- 
responding to a pushing beam intensity 1/4 of a MOT 
beam intensity), the flux increases sharply and then falls 
off with increasing power (see fig. 5). The simulation pre- 
dicts exactly the same critical power, without adjustable 
parameters (see fig. 5). The decrease at higher power can 
be understood if one examines the simulated divergence 
of the atomic jet, which grows with increasing pushing 
beam power. This effect is attributed to an insufficient 
short transverse cooling time due to the strong acceler- 
ation (see discussion below). Experimentally the jet ve- 
locity is deduced from measurements like flg. 2. With in- 
creasing pushing beam power it grows from 12 to 15m/s 
with an average width of lOm/s. In the simulation, we 
find a smaller width of 1 m/s. This discrepancy is proba- 
bly due to the fact that we have completely neglected the 
heating due to spontaneous emission. The longitudinal ve- 
locity width is larger than that of the LVIS or 2D MOT; 
however, for the purpose of loading a recapture MOT the 
velocity width does not matter. 

Detuning. The complex behaviour of the flux on the push- 
ing beam detuning (Jpush) is qualitatively very well repro- 
duced by the simulation (see figs. 6 and 7). If the pushing 
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Fig. 5. Dependence of the atomic flux on the pushing beam 
power. Flux is normalized, see text. The dots are experimental, 
the solid line is simulation. 
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Fig. 7. Simulation of the dependence of the atomic flux on 
pushing beam detuning. The flux is normalized as indicated in 
the text. 
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actual beam inclination and polarisation make the sit- 
uation a bit more complicated). For an atom traveling 
downwards in the extraction column, the |e,m = — 1) 
level approaches the MOT beam resonance at negative 
detuning. At the same time, the |e,m = +1) level ap- 
proaches the pushing beam resonance at positive detun- 
ing. When (5push — I^motIi the accelerating pushing beam 
and the decelerating MOT beams stay equally close to 
resonance throughout the extraction, and the atoms leave 
slowly. The extraction time is ^ 8 ms and the atoms are 
cooled transversely leading to a large recapture flux in the 
lower MOT. Finally if (5push > |<5mot|, the detuning of the 
|e,m = —1) level from the recentering MOT light is al- 
ways less than the detuning of the |e, m = +1) level from 
the pushing beam light, and so the trap is decentered but 
not destabilized (analogous to the behaviour at a large red 
detuning) . 



beam detuning is negative and exceeds the MOT beam 
detuning |(5push| > |<5mot|, the trap is decentered, but not 
yet leaking. Remember that the intensity of the pushing 
beam is about the same as for the MOT beams, so that 
as the detuning is increased the pushing radiation pres- 
sures becomes weaker than the trapping pressure. With 
zero or small blue detuning, atoms are resonantly acceler- 
ated, and their extraction is too fast to allow for efficient 
transverse cooling. These atoms leave at high divergence 
and are lost. Generally the simulation finds a 1 : 1 cor- 
relation of extraction time (flight time ~ 10 ms from the 
center of the trap to the depumping region) with diver- 
gence. Clearly, transverse cooling takes a certain time, and 
if the extraction acceleration is too strong, losses due to a 
high beam divergence are inevitable. 

For a blue detuning of the pushing beam such that 
^push — |i5mot|, a prominent peak in the flux appears in 
both the experiment and the simulation. To interpret this 
result we use the model of a \g, Jg = 0) ^ \e, Jg = 1) 
transition in a one dimensional magneto-optical trap (the 



Complementary Numerical Study: Waist. With a very 
small pushing beam size < 0.4 mm atoms drift out of the 
extraction column and are decelerated. They are recycled 
forever or leave the trap with high divergence. For a large 
spot size of > 1.5 mm atoms are not all extracted from the 
center and so many are not cooled sufficiently transversely. 
Both cases induce losses. 



5 Comparison and Conclusion 

Certainly, there are other techniques for the directed trans- 
fer of cold atoms from a VCMOT into a jet. A moving 
molasses launch Q provides a rather cold beam but low 
flux. A pulsed MOT launched by a resonant beampush is 
heated in the absence of transverse cooling beams ^ . Dur- 
ing the launch ~ vlOOO photons are spontaneously emit- 
ted into the transverse plane, while in continuous schemes 
there is transverse cooling during extraction. As a result 
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there is then no need for magnetic guiding , to achieve 
an elevated transfer efficiency. 

Continuous schemes suffer less from interparticlc in- 
teractions, since the steady state source cloud stays small. 
Leaking MOTs therefore accumulate atoms with the ini- 
tial capture rate of the MOT. The Atom Faucet provides 
a 50 % transfer efficiency from first capture, through the 
differential pumping tube, and to a recapture MOT in an 
UHV cell. It creates an extraction column that is typical of 
leaking MOT systems with a flexible design and without 
optical parts inside the vacuum chamber. 

The flux of = 1 X 10®atoms/s at a background 
vapour pressure of -PRbs? = 7.6 x 10~^mbar is equal to 
that of the low power version of the LVIS in Q and su- 
perior to the 20"*" MOT in this pressure region. The later 
design in turn provides very high flux at high pressure, 
since it minimizes the source trap density. We did not 
explore pressures that were incompatible with the UHV 
requirements in our recapture cell and found no devia- 
tion from the linear scaling of the flux with pressure up to 
pressures of 4 x 10^* mbar. Essentially, the Atom Faucet 
transplants to a MOT at 10^^^ mbar the loading rate of a 
MOT at few 10"^ mbar. 

We have also presented a 3D simulation of the atomic 
motion in multiple laser fields with an inhomogeneous 
magnetic field, neglecting interactions and fluctuations. 
We find that the transverse cooling inside the extraction 
column turns out to be a crucial element for the satis- 
factory performance of leaking MOT atom sources. Our 
simulation overestimates the capture rate, but predicts 
well the measured parameter dependences. Moreover, it is 
readily adapted to an arbitrary laser and _B-field configu- 
ration. 

We are indebted to F. Pereira dos Santos for coming 
up with the child's name and to the ENS Laser Cool- 
ing Group for helpful discussions. This work was partially 
supported by CNRS, College de France, DRET, DRED, 
and EC (TMR network ERB FMRX-CT96-0002). This 
material is based upon work supported by the North At- 
lantic Treaty Organisation under an NSF-NATO grant 
awarded to K.M. in 1999. W.W. gratefully acknowledges 
support by the Studienstiftung des deutschen Volkes and 
the DA AD. 

Note added: After this work was completed, we became 
aware that a similar setup has been successfully achieved 
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